The Bimodal Colors of Centaurs and Small Kuiper Belt Objects 
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Abstract. Ever since the very first photometric studies of Centaurs and Kuiper Belt Objects (KBOs) their 
visible color distribution has been controversial. That controversy gave rise to a prolific debate on the origin 
of the surface colors of these distant icy objects of the Solar System. Two different views attempt to interpret 
and explain the large variability of colors, hence surface composition. Are the colors mainly primordial and 
directly related to the formation region, or are they the result of surface evolution processes? To date, no 
mechanism has been found that successfully explains why Centaurs, which are escapees from the Kuiper 
Belt, exhibit two distinct color groups, whereas KBOs do not. In this letter, we readdress this issue using a 
carefully compiled set of B - R colors and Hr(o) magnitudes (as proxy for size) for 253 objects, including 
data for 10 new small objects. 

We find that the bimodal behavior seen among Centaurs is a size related phenomenon, common to both 
Centaurs and small KBOs, i.e. independent of dynamical classification. Further, we find that large KBOs 
also exhibit a bimodal behavior of surface colors, albeit distinct from the small objects and strongly 
dependent on the 'Haumea collisional family' objects. When plotted in B - R, //«(ff) space, the colors of 
Centaurs and KBOs display a peculiar N shape. 
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1. Introduction 



Discovered just 20 years ago ( Jewitt & Lu u|1993) , the Kuiper 
Belt holds a vast population of icy bodies orbiting the Sun be- 
yond Neptune. Stored at very low temperatures (~30-50 K), the 
Kuiper Belt Objects (KBOs) are expected to be well-preserved 
fossil remnants of the solar system formation. Presently, ~1600 
KBOs have been identified and classified into several dynam- 
ical families (see Appendix [A] and Gladman et al. 2008 for 
a review). KBOs which dynamically evolve to become Jupiter 
Family Comets (JFCs) form a transient population, the Centaurs, 
with short-lived chaotic orbits between Jupiter and Neptune 
( |Kowal et al.IT977)|Fernandez|1980HLevison & Duncan|1997| l. 

Between 1998 and 2003, we witnessed a debate on the sur- 
face colors of KBOs and Centaurs. One team used very accurate 
surface colors and detected that KBOs were separated into two 
distinct color groups ( |Tegler & Romanishin| 1998) [2000] [2003 1. 
Other teams did not find evidence for such color bimodality 
(Baru cci et aT][T999l |Jewitt & Luul|200Tj |Hainaut & Del santi 



2002| i. Careful reanalysis of the data by Peixinho et al. (2003) 



indicated that only the Centaurs display bimodal colors, i.e. they 
are distributed in two distinct color groups, one with neutral 
solar-like colors, and one with very red colors. KBOs on the 
other hand exhibit a broad continuous color distribution, from 



neutral to very red, with no statistical evidence for a color gap 
between the extrema (Te gler et al.|2008j for a review). 

The relevance of this controversy lays on two possible in- 
terpretations: i) KBOs and Centaurs are composed of intrinsi- 
cally different objects, with distinct compositions, which proba- 
bly formed at different locations of the protosolar disk, ii) KBOs 
and Centaurs are originally similar but evolutionary processes al- 
tered them differently, hence their color diversity. Most research 
focused on the latter hypothesis, offering little improvement on 
our understanding of the color distributions. Luu & Jewitt ( 1996) 
proposed that the competition between a reddening effect of 
irradiation of surface ices ^Thompso n et al.||1987| l and a blu- 
ing effect due to collisional induced resurfacing of fresh, non- 
irradiated, ices might generate the observed surface colors. The 
same authors, however, rejected this model as being the primary 
cause of the color diversity, due to the lack of predicted rota- 
tional color variations (Jewitt & Luu 2001). Based on the same 
processes, Gil-Hutton (2002) proposed a more complex treat- 
ment of the irradiation process, by implying an intricate structure 
of differently irradiated subsurface layers. However, the colli- 
sional resurfacing effects became very hard to model, thus mak- 
ing it very hard to provide testable predictions. Later, Thebault 
|& Doressoundiram] ( |2003| l showed that the collisional energies 
involved in different parts of the Kuiper Belt did not corrobo- 



1 



rate the possible link between surface colors and non-disruptive 
collisions. 



Table 1 : Filters specifications 



Delsanti et al. ( 2004 1 refined the first-mentioned model by 
considering the effects of a possible cometary activity trig- 
gered by collisions, and a size/gravity-dependent resurfacing. 
Cometary activity can modify the surface properties through 
the creation of a neutral-color dust mantle. |Jewitt ((2002 ) sug- 
gested that this process could explain why no JFCs are found 
with the ultra-red surfaces seen in about half of the Centaurs. 
It has also been proposed that the sublimation loss of surface 
ice from a mixture with red materials may be sufficient to make 
the red material undetectable in the visible wavelengths (Grundy 
|2009| . These might explain the Centaur color bimodality, as long 
as all were red when migrating inwards from the Kuiper Belt. 
Although promising, these models did not provide an explana- 
tion for the color bimodality of Centaurs, as they fail to repro- 
duce the bluest colors observed and their frequency. 

2. Motivation for This Work 

We find it striking that the objects with both perihelion 
and semi-major axis between Jupiter and Neptune's orbits, the 
Centaurs — by definition — , would display a different color dis- 
tribution than physically and chemically similar objects with a 
semi-major axis slightly beyond Neptune's orbit, as is the case 
for Scattered Disk Objects (SDOs), for instance, or any other 
KBOs. There is no evident physical consideration that would ex- 
plain the apparently sudden 'transition' in surface color behavior 
(from bimodal to unimodal) precisely at Neptune's orbital semi- 
major axis on =30.07 AU. This difference between Centaurs and 
KBOs is particularly puzzling since there is neither a sharp dy- 
namical separation between them, (the definition is somewhat 
arbitrary), nor a clearly identified family of KBOs in their origin. 
Although SDOs are frequently considered as the main source of 
Centaurs, Neptune Trojans, Plutinos, and Classical KBOs have 
been demonstrated as viable contributors (|Horner & Lykawka 
[20T0l |Yu & Tremainel[T999) |Volk & Malhotra||2008[ respect 
tivelly). Further, Centaurs possess short dynamical lifetimes of 
~ 5 ■ 10 5 -3 ■ 10 7 yr before being injected as JFCs or ejected again 
to the outer Solar System ( |Horner et al.|2004| . If some surface 
evolution mechanism, dependent on heliocentric distance, is re- 
sponsible for the bimodal behavior of Centaurs, it must be act- 
ing extremely fast such that no intermediate colors are ever seen 
among them. Besides surface color bimodality, the most distinc- 
tive characteristic of Centaurs compared to 'other' KBOs is their 
small size. Known KBOs are mostly larger than Centaurs, sim- 
ply because they are more distant and thus smaller objects are 
harder to detect. 

In this work, we address the issue of the color distributions 
of Centaurs and KBOs. We present new data on seven intrin- 
sically faint (thus small) KBOs and three Centaurs, combined 
with a new compilation of 253 published B-R colors, and avail- 
able m«(l, I, a) magnitudes, or Hr{o), i.e. absolute magnitude 
non-corrected from phase effects, and some identified spectral 
features. We study this large sample of colors (including ob- 
jects from all dynamical families) versus absolute magnitude as 
a proxy for size, with the implicit assumption that surface colors 
are independent of dynamical classification. We present the most 
relevant results, found in B - R vs. Hr{o) space. 



3. Observations and Data Reduction 

Observations of 7 KBOs and 1 Centaur were taken at the 
8.2 m Subaru telescope, on 2008-07-02, using 0."206/pix 





8.2m Subaru 


UH 2.2m 
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Center Width 


Center Width 


B 


4400 1080 


4480 1077 


R 


6600 1170 


6460 1245 



FOCAS camera in imaging mode with 2x2 binning (2 CCDs of 
2048 x 4096 pixels, Kashika wa"et al.|2002| l. Weather was clear 
with seeing ~0.7". We used the University of Hawaii UH 2.2 
m telescope, to observe 2 Centaurs on 2008-09-29, with the 
0."22/pixel Tektronix 2048 x 2048 pixels CCD camera. Weather 
was clear with seeing ~0.9". Both telescopes are on Mauna Kea, 
Hawaii, USA. Images from both instruments were processed 
using IRAF's CCDRED package following the standard tech- 
niques of median bias subtraction and median flat-fielding nor- 
malization. 

Standard calibration was made observing Landolt standard 
stars ( |Landolt|l 992 ) at different airmasses for each filter, obtain- 
ing the corresponding zeropoints, solving by non-linear least- 
square fits the transformation equations, directly in order of R 
and (B - R\ using IRAF's PHOTCAL package. The charac- 
teristics of the filters used on each telescope were essentially 
equal (Tab. [TJ. Subaru's data was calibrated using Landolt stan- 
dard stars: 107-612, PG1047+003B, 110-230, Mark A2, and 
1 13-337, taken repeatedly at different airmasses. UH2.2m's data 
was calibrated, analogously, using the stars: 92-410, 92-412, 
94-401, 94-394, PG2213-006A and PG2213-006B. These stars 
have high photometric accuracy and colors close to those of the 
Sun. We have used the typical extinction values for Mauna Kea, 
k B = 0.19, and k R = 0.09 ( |Krisciunas et~aT1[T987l and CFHT 
Info Bulletin #19). All fits had residuals rms < 0.02, which were 
added quadratically to the photometric error on each measure- 
ment. Targets were observed twice in B and twice in R bands, 
to avoid object trailing in one long exposure. Each two B or R 
exposures were co-added centered in the object, and also co- 
added centered on the background stars. The former were used 
to measure the object, the latter to compute the growth-curve 
correction. The time and airmass of observation were computed 
to the center of the total exposure. We applied growth-curve 
correction techniques to measure the target's magnitudes using 
IRAF's MKAPFILE task (for details, see |Peixinho et al.|200 4). 
Observation circumstances and results are shown in Table [2] 



4. Compilation of Data 

We compiled the visible colors for 290 objects (KBOs, 
Centaurs, and Neptune Trojans) for which the absolute magni- 
tude in R or V band was accessible {e.g. with individual mag- 
nitudes and observing date available), and surface spectra in- 
formation for 48 objects, as published in the literature to date 
(Feb. 2012). We computed the absolute magnitude H R (a) = 
mail, I, a) = R-5 log(r-A), where R is the R-band magnitude, r 
and A are the helio- and geocentric distances in AU, respectively. 
In this compilation, 253 objects have B-R color available which 

and 48 have also spec- 



is the focus of this paper (see Table A. 1 
tral information. The description of the compilation method is 
presented in Appendix [A] Sun-Object-Earth phase angles a are, 
typically, less than 1.5° for KBOs and less than 4° for Centaurs. 
Measurements of magnitude dependences on the phase angle for 



2 



Table 2: Observational circumstances and photometric results of this work's data 



Object 


Dyn. Class* 


Telescope 


UT Date 


r[AU] 


A[AU] 


a[°] 


R 




B-R 


H R (a) 


(130391) 


2000 JG 8 i 


2:1 


Subaru 


20080702UT07:24:58 


34.073 


34.817 


1.2 


23.12±0.03 


1.42±0.06 


7.75+0.06 


(136120) 


2003 LG 7 


3:1 


Subaru 


20080702UT09:42:53 


32.815 


33.659 


1.0 


23.544 


-0.05 


1.27+0.09 


8.32+0.05 


(149560) 


2003 QZ 91 


SDO 


Subaru 


20080702UT13:08:33 


25.849 


26.509 


1.7 


22.484 


:0.03 


1.30+0.05 


8.30+0.03 




2006 RJ IOT 


Nep. Trojan 


Subaru 


20080702UT14:07:50 


30.760 


30.534 


1.9 


22.274 


.0.02 


1.90+0.04 


7.40+0.02 




2006 SQ372 


SDO 


Subaru 


20080702UT1 1:45:34 


23.650 


24.287 


1.9 


21.554 


0.02 


1.78+0.03 


7.71+0.05 




2007 JK43 


SDO 


Subaru 


20080702UT08:08:13 


23.113 


23.766 


1.9 


20.734 


:0.02 


1.40+0.03 


7.03+0.02 




2007 NC 7 


SDO 


Subaru 


20080702UT1 1:30:49 


20.090 


20.916 


1.7 


21.194 


:0.02 


1.28+0.03 


8.07+0.02 


(281371) 


2008 FC 76 


Cent 


Subaru 


20080702UT11:13:05 


11.119 


11.793 


3.8 


19.79 = 


t0.02 


1.76+0.02 


9.18+0.04 




2007 RH 28 i 


Cent 


UH2.2m 


20080929UT12:43:47 


17.081 


17.956 


1.6 


20.85 = 


t0.03 


1.20+0.05 






2007 RH 283 


Cent 


UH2.2m 


20080929UT12:57:51 


17.081 


17.956 


1.6 


20.90+0.03 


1.28+0.06 






mean... 


















1.24+0.07 


8.44+ 0.04 




2007 UM 126 


Cent 


UH2.2m 


20080929UT08:56:52 


10.191 


11.177 


0.9 


20.434 


:0.03 


1.21+0.05 






2007 UM 126 


Cent 


UH2.2m 


20080929UT09:06:41 


10.191 


11.177 


0.9 


20.534 


.0.03 


0.92+0.04 






2007 UM 126 


Cent 


UH2.2m 


20080929UT09:16:17 


10.191 


11.177 


0.9 


20.384 


:0.02 


1.12+0.04 






mean... 


















1.08+0.10 


10.16+0.04 



Notes. * Dynamical classes are: Centaur, Scattered Disk Object (SDO), Neptune Trojan (object located in 1:1 mean motion resonance with 
Neptune), 2:1, and 3:1, (objects located in 2:1 or 3:1 mean motion resonance with Neptune, respectively). For details on the classification see 
Appendix [A| 



these objects, i.e. phase coefficients y6[mag/°], are scarce but, so 
far, do not show evidence for extreme variability presenting an 
average value of/? = 0.11 + 0.05 (Belskaya et al. 2008). From 
the linear approximation Hr(u = 0°) » Hr(<x) — aB, by not cor- 
recting the absolute magnitude from phase effects we are slightly 
overestimating it. We will deal with this issue in Sec. [5] 

Recent works have shown that there is no strong corre- 
lation between object diameter D and geometric albedo pv, 
nor between geometric albedo pv and absolute magnitude Hr 
(|Stansberry et al.|20 08; Sant os-Sanz et al.|2012 Vilenius et al. 
|2012| |Mommert et al.||2012| i. However, from the 74 diameter 
and albedo measurements of Centaurs and KBOs made using 
Herschel and/or Spitzer telescopes, published in the aforemen- 
tioned works, we verify that Hr and D correlate very strongly 
with a Spearman-rank correlation of p = -0.92^2, with 
a significance level SL <k 0.01% (error bars computed us- 
ing bootstraps, for details see |Doressoundiram et al. 2007 1. 
Consequently, absolute magnitude is a very good proxy for size. 



5. An TV-shaped Doubly Bimodal Structure 

In Fig. [T] we plot /?-band absolute magnitude Hr(u) (proxy 
for object's size) against B-R color for all (n = 253) objects in 
our database. The cloud of points forms a recognizable N shape 
with an apparent double bimodal structure in color. The smaller 
objects (upper part of the plot) show a bimodal B-R distribu- 
tion. Although apparently dominated by Centaurs, this bimodal 
distribution also includes KBOs of similar Hr(o), which sug- 
gests that the bimodal structure in B-R color is a property of 
the smaller objects in general, regardless of their dynamical 
family. This bimodality appears do disappear for objects with 
Hr{u) < 7 where the B-R color distribution seems unimodal. 
Most interestingly, we note that towards the larger objects (lower 
part of the plot) the colors suggest the presence of another bi- 
modal behavior, with the gap between the two groups shifted 
towards the blue with respect to the 'small' object bimodality. 
This new Targe' object bimodality is explicitly reported for the 
first time. 

When performing hypotheses testing one should adopt a crit- 
ical value of significance a. The value a is the maximum proba- 
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Fig. 1: B - R vs. H R (a) plot of all 253 objects. KBOs are represented 
by solid circles and Centaurs by white dotted solid circles. Objects with 
Hn( a ) ^ 6.8 separate into two color groups with a 'gap' centered at 
B — R ~ 1.60. Objects with H R (a) ^ 5.0 also show statistical evidence 
for separation in two colors groups but with a 'gap' centered at B — 
R ~ 1.25. Objects spectra with known features of water ice, methane, 
methanol, and featureless spectra, are coded using colors as described 
in the legend. There is no obvious/clear connection between B-R colors 
and the presence of spectral features. 



bility (risk) we are willing to take in rejecting the null hypothesis 
Ho {i.e. to claim no evidence for bimodality) when it is actu- 
ally true (i.e. data is truly bimodal/multimodal) — also called 
type I error probability. Such value is often a source of debate, 
as the theories of hypotheses testing themselves (e.g. Lehmann 
|1993| ). The decision relies mostly whether the effects of a right or 
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wrong decision are of practical importance or consequence. The 
paradigm is: by diminishing the probability of wrongly reject a 
null hypothesis (e.g. decide for bimodality when bimodality was 
not present in the parent population) we increase the probabil- 
ity of wrongly accepting the null hypothesis (i.e. deciding for 
unimodality when bimodality was in fact present), also called 
type II error probability, or risk factor B. Some authors and/or 
research fields, consider that there is only sufficient evidence 
against Hq when the achieved significance level is SL < 0.3%, 
i.e. using a = 0.3% (the 3<x Gaussian probability), others require 
even a = 0.0003% (6<x). Such might be a criterion for rejec- 
tion of Hq but not a very useful 'rating' for the evidence against 
Hq, which is what we are implicitly doing. We rate the evidence 
against Hq following a most common procedure in Statistics: 
SL < 5% — reasonably strong evidence against Hq, S L < 2.5% 
— strong evidence against Hq, and SL < 1% — very strong ev- 
idence against Hq (e.g. Efron & Tibshirani |1993| l, adding also 
the common procedure in Physics: SL < 0.3% — clear evidence 
against Hq. Further, for better readability, throughout this work 
we may employ the abuse of language 'evidence for bimodal- 
ity' instead of the statistically correct term 'evidence against uni- 
modality'. 

Using the R software's (versi on 2.14.1; |R Deve lopment 
Core Team|201 1\ Dip Test pack age ( |Hartigan|19 85 ; Hartiganj&j 



Hartigan 1985; Maechler 20 1Tb we test the null hypothesis Hq: 
'the sample is consistent with an unimodal parent distribution' 
over all objects in the B - R vs. H R (a) space, against the alter- 
native hypothesis H\. 'the sample is not consistent with an uni- 
modal parent distribution' (hence it is bimodal or multimodal). 
The full sample, in spite of the apparent two spikes, shows no rel- 
evant evidence against color unimodality, neither with (n = 253, 
SL = 17%) nor without (n = 224, SL = 41%) Centaurs (see 
Fig. [2] a). The Centaur population (n = 29) shows strong ev- 
idence against unimodality at 1.6%. Removing the 3 brightest 
Centaurs (with H R (a) > 6.6) improves the significance to 0.3%. 
To refine the analysis and test different ranges in H R (a) we ran 
the Dip Test on sub-samples using a running cutoff in H R (a) that 
was shifted by 0.1 mag between consecutive tests. 



Bimodal distribution of 'small' objects: We performed iter- 
ative Dip Tests with a H R:cut starting at the maximum H R (a) 
value, and decreasing in steps of 0.1 mag; in each iteration we 
run the test on those objects above the cutoff line (i.e. with 
H R (a) ^ H Rcut ). We stop shifting H R:cu , when we detect the 
maximum of evidence against unimodality (i.e. a minimum of 
significance level, henceforth accepting the alternate hypothesis 
'the distribution is bimodal/multimodal' ) Evidence for bimodal- 
ity at significance levels better than 5% start to be seen for ob- 
jects with H R (a) ^7.1. This evidence peaks at a significance of 
0. 1 % for the 124 faint objects with H R (a) > 6.8. 

We propose that the visible surface color distribution of 
(non-active) icy bodies of the outer Solar System depends only 
on objects size, and is independent of their dynamical classifi- 
cation. No mechanism has yet been found to explain the color 
bimodality only for Centaurs. However, since such mechanism 
might exist even if not yet found, we re-analyze the sample re- 
moving the Centaurs. Naturally, the sampling of the smaller ob- 
jects diminishes considerably, hence reducing the statistical sig- 
nificance against the null hypothesis (i.e. increases the probabil- 
ity of observing two groups on a purely random distribution of 
colors). Nonetheless, the 98 remaining objects with H R (a) ^ 6.8 
show evidence for bimodality at a significance level of 3.5%, 
reaching a significance minimum of 1.8% for the 165 objects 



with H R (a) 5.8. In both cases the 'gap' is centered around 
B -R ~ 1.60 (see Figs. [T]and|2]b). 

Bimodal distribution of 'large' objects: We test the bright- 
est part of the sample using a cutoff limit starting at the min- 
imum H R (a) value; we consider objects below the cutoff (i.e. 
brighter than H R:cut ) and shift it up in steps of 0.1 mag. We 
find very strong evidence against unimodality for objects with 
H R (a) < 5.0 (SL = 0.9%). Data still shows reasonably strong 
evidence against unimodality for objects up to H R (a) < 5.6. 
The 'gap' is located at B - R ~ 1.25. There are no Centaurs 
in this brightness range. Explicitly, evidence for 'large' objects 
bimodality has not been previously reported, (see Figs. [T] and 
[2] c). Removing from the sample the 7 objects belonging to 
the 'Haumea collisional family' ( Brown et al.|2007b} |Snodgrass 
|et al.|2010) , all clustered on the lower left 'leg' of the N shape, 
erases the statistical evidence against the null hypothesis, even if 
still suggestive to the eye. Therefore, with the present data sam- 
ple, the 'evidence for bimodality' among bright KBOs cannot be 
stated as independent from the peculiar properties of the Haumea 
collisional family. 

The 'intermediate' size continuum: The 91 objects with 6.8 > 
H R (a) > 5.0, which include 3 Centaurs, do not show evidence 
against a unimodal behavior (SL - 98.0%) even if a small gap 
seems suggestive to the eye (see Figs.[T]and|2]d). However, sta- 
tistically, their inclusion in the fainter group does not decrease 
the significance below the 'strong evidence against unimodal- 
ity', i.e. SL = 2.5% (see Figs. [T] and [2] d). On the other hand, if 
added to the 'large' objects the statistical evidence for bimodal- 
ity of 'large' objects does not hold. 

To check for the effects of non-correcting H R (a) from phase 
angle effects we performed Monte-Carlo simulations. First, we 
compute all the possible a values and their probability distri- 
bution for an 'average' Centaur with semi-major axis a — 15 
AU. The maximum a is 3.8° being the median value 3.2°. 
Analogously, we do the same for a KBO with a — 40 AU. 
The maximum a is 1.4° and the median value 1.2°. Therefore, 
on average, our absolute magnitudes might be overestimated by 
AH R ~ 0.35, for Centaurs, and by AH R « 0.13, for KBOs. 
Simulating 1000 'phase-corrected' H R data-samples, following 
the probability distribution of the corresponding a angles did 
not alter any of the results obtained using simply H R (a). 



6. Interpretation 

Our analysis shows that the B - R colors of Centaurs and 
KBOs when plotted as a function of H R (a) display an N-shaped, 
double bimodal behavior. The color distribution seems to de- 
pend on object size (intrinsic brightness) instead of dynam- 
ical family. Using the brightness-size-albedo relation D^,, = 
2 V2.24 ■ 10 16 ■ W )A( - H «°- H ^ I p R , with solar H RB = -27.10, the 
main issue is to choose a canonical geometric albedo value p R . 
Recent works (Stansber ry et al.||2 008, Santos-S anz et al.||2012 



|Vilenius et al.|2012||Mo mmert e t al.| 2012| l show a wide range 



of values, for each dynamical family, in some cases far from 
the 0.04 value previously assumed from comet studies. As we 
need only a rough estimate of size ranges, we pick the average 
value of p R — 0.09. Using this parameter, objects with diame- 
ters 165 < Dkm S 380 present a rather continuous range of B-R 
colors. 
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Fig. 2: Histograms of B — R colors from selected H R (a) ranges: a) All 
the 253 objects. Taken globally do not exhibit statistical evidence for bi- 
modality, which was known to exist among Centaurs, b) The 124 'small' 
objects, with H R (a) 6.8. Evidence for bimodal behavior is clear and 
still present when removing Centaurs, c) The 38 'large' objects, with 
Hr{o) ^ 5.0. A bimodal behavior is shown but it loses the statistical 
significance without the 'Haumea collisional family' objects, d) The 91 
'intermediate' size objects, 6.8 > H R {a) > 5.0. Regardless of the ap- 
parent small gap at B - R ~ 1.3 there is no statistical evidence for two 
separate groups. 



Visible and near-infrared (NIR) spectroscopy for about 75 
bright objects ( Barucci et al.|20l"T| for a review) also indicates 
that the surface compositions of KBOs and Centaurs is very di- 
verse. The largest objects are coated in methane ice, while inter- 
mediate size objects display water-ice features, sometimes with 
traces of other volatiles. Small KBOs generally have featureless 
spectra. The presence of volatiles on the surface of an object 
may be related to its ability to retain them, i.e. to its size and 
temperature (Schaller & Bro wn|2007|l. It should also depend on 
the subsequent irradiation history ( Brown et al.|201 1 1. However 
no correlation can be made to date between visible colors and 
NIR spectral properties. For example, two objects of compara- 
ble size, Quaoar and Orcus, both exhibit water ice-dominated 
surfaces but have, respectively, very red and neutral visible col- 
ors ( |Delsantietal.|2010[ ). 

Objects smaller than ~ 100- 150 km, including most of the 
known Centaurs, are believed to be fragments from the collision 
of larger objects (|Pan & Sari|2005[). Predicting the properties of 



these fragments is a complex task, but the field shows promising 
advances (for a review, see |Leinhardt et al.|2008| >. An immediate 
hypothesis is that red and neutral objects are the only possible 
outcomes of a disruptive collision. Thermal evolution modeling 
suggests that KBOs, especially large ones, should have a lay- 
ered structure, including some liquid water leading to a complete 



differentiation of the object (Merk & Prialnik 2006 Guilbert- 
|Lepoutre et al.|20lT) . A catastrophic collision could result in the 
formation of fragments with very different properties, depending 
on whether they come from the core of the parent body, or its 
mantle, or some subsurface layers. However, our current knowl- 



edge of KBOs internal properties and evolution is still incipient 
to support or discard such an hypothesis. Besides, it is hard to 
understand why objects with B-R ~ 1.6 (in the gap of the 
small object's bimodal distribution) should be inexistent. Maybe 
their relative number is so small compared to the neutral and 
red groups that we can hardly observe them, leading to another 
puzzling question. Research on these aspects should be encour- 
aged, in particular the detection and measurement of many more 
small objects — KBOs and Centaurs — could help further con- 
straining their color distribution and other properties. The ob- 
jects in the 'intermediate' H R (a) range (6.8 > H R (a) > 5.0) 
seem unimodally distributed in B - R color; they might repre- 
sent a transition phase between the two bimodal distributions. 
These medium-sized objects are probably too large to be rem- 
nants from disruptive collisions, and too small to have recently 
undergone cryovolcanic activity (they may not even be differen- 
tiated). They might, actually, represent the only group where the 
outcomes of the combined effects of different birthplaces, space 
weathering and thermal processing can be studied or analyzed. 

The evidence for bimodal distribution among the largest ob- 
jects is also puzzling. These are supposedly the best studied ob- 
jects, yet the evidence for a bimodal distribution of their surface 
colors has never been reported. Nonetheless, when removing the 
7 Haumea collisional family objects from our sample it no longer 
provides evidence against an unimodal distribution, even if ap- 
parent to the eye. This issue should be further analyzed when 
larger sampling is available. 

In this work, we confirm that there is no noticeable link be- 
tween the surface composition of an object and its visible col- 
ors. Objects hosting water ice are distributed both among large 
and small objects, and among red and blue ones. When it comes 
to volatiles such as methane (CH4) or methanol (CH3OH), we 
find that they are also distributed among all groups, although 
they might be more difficult to detect on small/fainter objects. 
We nonetheless find a cluster of featureless objects among the 
red group of large objects: these might represent the most ir- 
radiated/oldest surfaces in the overall population. Therefore, it 
seems that a simple explanation such as the model of atmo- 
spheric escape proposed by Schaller & Brown ( 2007 ) might not 
be sufficient to explain the colors and compositions of KBOs. 
The reason why they evolved in two different color groups can be 
very complex, and should involve different thermal, collisional, 
irradiation histories, on top of possible different birthplaces. 

7. Summary 

In this work we analyze the B-R color distribution as a func- 
tion of H R (a) magnitude for 253 Centaurs and KBOs, including 
10 new measurements, and with the information on their NIR 
spectral features. Using the known diameters, D, and albedos, 
Pv, of 74 of these objects we verify that H R and D correlate 
very strongly (p = -0.92+° ^, SL <sc 0.01%) validating H R as a 
good proxy for size. Further, through simulations, we show that 
not correcting H R (a) to H R (a = 0°) does not change any of the 
global results. Our analysis shows: 

1. The B-R vs. H R (a) color distribution is Af-shaped, evidenc- 
ing that B-R colors are probably dominated by a size effect 
independent from dynamical classification. 

2. Small objects, including both KBOs and Centaurs, display 
a bimodal structure of B - R colors at 0.1% significance 
level (i.e. objects with H R (a) ^ 6.8, or < 165, assum- 
ing p R = 0.09) with the 'gap' centered at B - R ~ 1.60. 
Removing Centaurs from the sample reduces greatly the 



5 



sampling on small objects reducing also the significance of 
the result to 3.8%. 

3. Large objects evidence also for a bimodal structure, with 
minimum significance of 0.9%, for Hr(o) < 5.0 (D^m > 
380, assuming = 0.09), and color 'gap' centered at 
B - R ~ 1.25. Reasonable evidence for this bimodality 
starts when considering only objects with Hr(o) < 5.6 
(Dkm <; 290) dropping below the critical 5% when reaching 
Hr(c?) < 4.4 (Dk m ^ 500). However, this behavior seems 
dominated by the presence of 7 Haumea collisional family 
objects which 'cluster' at the lower left leg of the /V-shape. 
Once removed, there is no statistical evidence against com- 
patibility with a random unimodal distribution for the larger 
KBOs. 

4. Intermediate sized objects do not show incompatibility with 
a continuum of B - R colors (i.e. 6.8 > H R (a) > 5.0, or 
165 < Dtm tS 380, assuming /?« = 0.09). These objects 
seem too large to be remnants from disruptive collisions and 
too small to hold cryovolcanic activity. They might be the 
best targets to study the combined effects of different birth- 
places, different space weathering, and different thermal pro- 
cessing. Further studies are encouraged. 

5. Inspecting the NIR spectral properties against B - R colors 
shows no obvious link between the colors and the chemical 
composition of the objects' surfaces. 
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Appendix A: Compiled Database 

For each object, we compute the average color index from 
the different papers from data obtained simultaneously in B and 
R bands {e.g. contiguous observations within a same night). 
When individual R apparent magnitude and date is available, we 
compute the H^ia) — R - 5 log (r ■ A), where R is the R-band 
magnitude, r and A are the helio- and geocentric distances at 
the time of observation in AU, respectively. When V and V — R 
color is available, we derive an R and then Hr{o) value. We do 
not correct for the phase angle a effect as we need only a gen- 
eral estimation of the absolute magnitude for our complete sam- 
ple. In addition, few objects have phase correction coefficient 
available in the literature, and no universally accepted canoni- 
cal values per dynamical class can be strictly adopted. Table A. 1 



presents the resulting values. This table includes also spectral 
information on the presence of water ice, methanol, methane, or 
confirmed featureless spectra, as available in the literature. We 
highlight only the cases with clear bands on the spectrum which 
were reported/confirmed by some other work. 

There is no strict definition for the dynamical classes of 
Centaurs and KBOs. Roughly speaking: objects orbiting in mean 
motion resonances with Neptune are called 'resonants' (if lo- 
cated in the 1:1 resonance are also known as Neptune Trojans, 
and known as Plutinos if located in the 3:2 resonance); Centaurs 
are the objects with orbits between those of Jupiter and Neptune; 
Scattered Disk Objects (SDOs), are those within probable grav- 
itational influence of Neptune; Detached KBOs, are those be- 
yond past or future gravitational influence by Neptune; Classical 
KBOs , are those with rather circular orbits beyond Neptune and 
below the 2:1 resonance region (being called Hot if their orbital 
inclination is higher than 5° or Cold if lower). 

To determine the dynamical class we first gathered the orbital 
elements, with epoch 201 1-12-05, from 'The Asteroid Orbital 
Elements Database', astorb.datpl maintained by the 'Lowell 
Observatory' based on astrometric observations by the 'Minor 
Planet Center'. Then, using the particular classification scheme 
suggested by Lykawka & Mukai (2007), including their analysis 
of objects located in the mean motion resonances (MMR) with 
Neptune, dynamical class was determined following all steps 
algorithm: 



10. i < 5° A { [q ^ 37 AU A (37 AU s? a ^ 40 AU)] V [q ^ 
3SAU A (42 AU s$ a 48 AU)] } => Cold Classical KBO 
(cCKBO) 

11. i ^ 5°Aq ^ 31AUA(31AU < a < 48 AU) => Hot Classical 
KBO (hCKBO) 

being q and a the object's perihelion and semi-major axis, re- 
spectively. Jupiter semi-major axis is a j, and Neptune's is a^. 
Note that throughout the algorithm an object can be reclassified. 

We are aware that there are more complex classification 
schemes, which may be more refined, but the boundaries be- 
tween families do not change significantly. We chose this one 
for its computational simplicity. 

References for the colors presented in Table A.l are : (1) 
I Luu & Jewittj (fl996| ); ( 2) |Lazzaro et al.| ([1997]) ; (3) 



|Martin et al. (2003); (4) |Romanishin et a 



1997); (5) 



Mart in et al.|([2002|l; (6)|Tegler & Romanishin|(|1998| l; (7)|Jewitt 



IRomanishin 



pOOO) ; (10 )|Delsanti et al.|(|2001) ; (1 1) 



Romanishin 

H1999); (14 ) |Boehnhardt eTaL] pM 
et al. (2007); (16) Doressoundiram et al 



Romon- 



Romon- 



& Luu| POOIFW [ Doressoundira met all (|2002| l; (9) 



Tegler & 



Tegler & 



( p997); (12)|Jewitt & Luu|(|1998| ); ( 13)|Barucci et al. 



(15) Doressoundiram 
(EUTJT); !(17)|Green etal 



( M997t ; ( 1 8) |Boehnhardt et al.|(|2002| >; (1 9)|Tegler & Romanishin 
pO Ol); (20) |Hainaut et al |(|2000|l; (21) |Sheppardj ( |2010| l; (22) 



Barucci et al. (2000); (23) Rabinowitz et al. (2008 



(24) Tegler 
l/survey.htm; 



et al. http)//www.physics. nau.edu/ tegler/re search 
(25)|Tegler et al.| ( [2003]l; (26)|Peixinho et"alJ(T200"TT >; (27)|Trujillo 

"")4); (29) |Ferrin et al. 

b); (31)|Santos-Sanz 
Potto etal IpOOl) ; (33) |Fornasier et al M 2004 
ram et al.| j|2005a| ); (35) |Gulbis e t al.| (|2006 



d2001|>; (30 ) IDoressoundiram et al. | (|2005bl 
etal.|H2009D;(32T g ' 
(34) Doressoundiram 



(36) |Rabinowitz etaL] d2007); (37) 



ISnodgrass et al. 



(2010); 



% 

(38) |Rabinowitz et al.| ( |2006); (39)|Lacerda et al.|(|2008) ; (40) 



Romanishin et al.l (12010 



(41) Sheppard & Trujillo (2006) 



References for the spectral features indicated in Table |A.l 
are : (a)|Rombn-Marfin et al.|<|2003]>; (b)|Cruikshank e"taLl(fl 998 V, 
(c)|Kern et al.|(|2000|); (d)|Guilbert et al.| (|2009b); (e)|Jewitt & 



|Luu|([2()()T); (f)|Brown et al.| ( |T999); (g) Barkum e et al. 
(h)|Guilbert et al | (|2009a|>; (i)IBarucci et al.|(|201 1|>; ?jj 
let al.|H2010); (k)|Grundy et al.|(|2005|l; (l)|Barucci et al. 



(2008) 



DeMeo 
q2010| ; 



(m) IDelsanti et al.| (|201 01)7(5) |Pinilla-Alonso et al.| <|2009); (o) 



[ Merlin et al.| ( |2009) ; (p) [Brown et al.| | |2007a) ; (q) |Schaller &| 
|Brown| ( |2008| ). 



1 . q < a j => Not analysed 

2. in 1 : 1 MMR with Neptune => Neptune Trojan 

3. in 3 : 2 MMR with Neptune => Plutino 

4. in other MMR with Neptune => Other Resonant 

5. q > aj A a < ^> Centaur 

6. aj < q < a N A a ^ a N => Scattered Disk Object (SDO) 

7. a N < q ^ 37 AU => Scattered Disk Object (SDO) 

8. O 40 At/Aa^ 48 At/=> Detached KBO (DKBO) 

9. 31 AU sC q < 40 AU Scattered or Detached KBO 
(SDKBO) 



ftp://ftp.lowell.edu/pub/elgb/astorb.dat.gz 
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Table A.l. Compilation of absolute magnitude H R (a), B - R colors, and spectral features used in this work 



Object 




Dynamical Class 


H R (a) 


B-R 


Spectral features References 


(2060) Chiron 


Centaur 


6.287±0.022 


1.010+0.044 


Water ice 


1, 2, 3, a 


(5145) Pholus 


Centaur 


7.198+0.056 


1.970+0.108 


Methanol 


4,b 


(7066) Nessus 


Centaur 


9.020±0.068 


1.847+0.165 




1 


(8405) Asbolus 


Centaur 


9.257+0.120 


1.228+0.057 


Water ice 


4, 5, c 


(10199) 


Chariklo 


Centaur 


6.569+0.015 


1.299+0.065 


Water ice 


6, 7, d 


(10370) 


Hylonome 


Centaur 


9.250+0.131 


1.153+0.081 




1,6,8 


(15760) 


1992 QB] 


Cold Classical 


6.867+0.121 


1.670+0.145 




1,7,9 


(15788) 


1993 SB 


Plutino 


8.032+0.122 


1.276+0.100 




7, 9, 10 


(15789) 


1993 SC 


Plutino 


6.722+0.074 


1.720+0.140 


Methane 


1,7, 11, 12, e 


(15810) 


1994 JRj 


Plutino 


6.867+0.077 


1.610+0.216 




13 


(15820) 


1994 TB 


Plutino 


7.527+0.091 


1.759+0.155 




1,7, 10, 11, 13 


(15874) 


1996 TL 66 


Scattered Disk Object 


5.131+0.144 


1.113+0.070 




6,7, 12, 13, 14 


(15875) 


1996 TP 66 


Plutino 


6.953+0.071 


1.678+0.123 




6, 7, 12, 13, 14, 15 


(15883) 


1997 CR 29 


Scattered Disk Object 


7.076+0.135 


1.260+0.128 




7, 16 


(16684) 


1994 JQ! 


Cold Classical 


6.618+0.117 


1.738+0.120 




17, 18, 19 


(19255) 


1994 VK 8 


Cold Classical 


7.016+0.163 


1.680+0.067 




9 


(19299) 


1996 SZ 4 


Plutino 


8.184+0.159 


1.299+0.102 




7, 9, 18 


(19308) 


1996 TO 66 


Resonant (19:11) 


4.530+0.044 


1.056+0.210 


Water ice 


6,7, 12, 13, 14, 20, 21, f 


(19521) 


Chaos 


Hot Classical 


4.442+0.069 


1.558+0.062 




8, 9, 10, 22 


(20000) 


Varuna 


Hot Classical 


3.345+0.059 


1.530+0.036 


Water ice 


8, g 


(20108) 


1995 QZ 9 


Plutino 


7.889+0.399 


1.400+0.050 




9, This work 


(24835) 


1995 SM 55 


Hot Classical 


4.352+0.040 


1.018+0.052 


Water ice 


8, 10, 14, 23, g 


(24952) 


1997 QJ 4 


Plutino 


7.389+0.114 


1.104+0.104 




7, 18 


(24978) 


1998 HJ 15 i 


Cold Classical 


7.008+0.050 


1.820+0.042 




19 


(26181) 


1996 GQ 21 


Resonant (1 1:2) 


4.467+0.090 


1.693+0.079 


Methanol 


18, 24, g 


(26308) 


1998 SM 165 


Resonant (2:1) 


5.757±0.119 


1.620+0.105 




9, 10, 15 


(26375) 


1999 DE 9 


Resonant (5:2) 


4.810+0.046 


1.536+0.056 


Featureless 


7, 8, 25, h 


(28978) 


Ixion 


Scattered Disk Object 


3.366+0.038 


1.634+0.035 


Water ice 


8, h 


(29981) 


1999 TDio 


Scattered Disk Object 


8.698+0.038 


1.230+0.028 


Water ice 


8, g 


(31824) 


Elatus 


Centaur 


10.439+0.107 


1.672+0.071 


Water ice 


8, 10, 26, g 


(32532) 


Thereus 


Centaur 


9.365+0.038 


1.190+0.032 


Water ice 


25, h 


(32929) 


1995 QY 9 


Plutino 


7.489+0.126 


1.160+0.150 




1, 13 


(33001) 


1997 CU 29 


Cold Classical 


6.173+0.078 


1.804+0.115 




7, 16, 22, 27 


(33128) 


1998 BU 48 


Scattered Disk Object 


6.889+0.127 


1.692+0.089 




8, 10 


(33340) 


1998 VG 44 


Plutino 


6.292+0.077 


1.511+0.055 




8, 14, 16, 24 


(35671) 


1998 SN 165 


Scattered Disk Object 


5.431+0.068 


1.123+0.082 




7, 10, 16 


(38083) 


Rhadamanthus 


Scattered Disk Object 


7.432+0.063 


1.177+0.109 




18 


(38084) 


1999 HB 12 


Resonant (5:2) 


6.718+0.050 


1.409+0.049 




16, 25, 27, 28 


(38628) 


Huya 


Plutino 


4.674+0.099 


1.539+0.062 


Featureless 


29, 7, 16, 18, g 


(40314) 


1999 KR I6 


Scattered Disk Object 


5.527+0.039 


1.872+0.068 




7, 18,27 


(42301) 


2001 UR I63 


Resonant (9:4) 


3.812+0.109 


2.190+0.130 


Featureless 


15, 30, 31, g 


(42355) 


Typhon 


Scattered Disk Object 


7.358+0.076 


1.292+0.071 


Water ice 


25, 28, h 


(44594) 


1999 OX 3 


Scattered Disk Object 


6.835+0.078 


1.839+0.087 


Water ice 


8, 9, 10, 15, 21,30, i 


(47171) 


1999 TC 36 


Plutino 


4.851+0.054 


1.740+0.049 


Water ice 


10, 16, 25, 32, h 


(47932) 


2000 GN171 


Plutino 


5.666+0.090 


1.559+0.066 


Featureless 


18, 24, h 


(48639) 


1995 TL 8 


Detached KBO 


4.667+0.091 


1.693+0.217 




8, 10,21 


(49036) 


Pelion 


Centaur 


10.157+0.112 


1.248+0.096 




9, 18 


(50000) 


Quaoar 


Hot Classical 


2.220+0.029 


1.588+0.021 


Methane 


25, 33, h 


(52747) 


1998 HMi 51 


Cold Classical 


7.417±0.100 


1.550+0.103 




19 


(52872) 


Okyrhoe 


Centaur 


10.775+0.078 


1.237+0.086 


Water ice 


10, 16, 32, g 


(52975) 


Cyllarus 


Centaur 


8.634+0.101 


1.803+0.102 




8, 10, 14, 25 


(53311) 


Deucalion 


Cold Classical 


6.662+0.060 


2.030+0.160 




27 


(54598) 


Bienor 


Centaur 


7.727+0.077 


1.158+0.075 


Methanol 


8, 10, 15, h 


(55565) 


2002 AWi 97 


Hot Classical 


3.156+0.059 


1.498+0.044 


Featureless 


24, 33, 34, h 


(55576) 


Amycus 


Centaur 


7.789+0.042 


1.814+0.044 


Water ic 


24, 28, 33, 34, i 


(55636) 


2002 TX 300 


Hot Classical 


3.296+0.047 


1.010+0.028 


Water ice 


25, 30, q 


(55637) 


2002 UX 25 


Scattered Disk Object 


3.486+0.084 


1.502+0.052 


Water ice 


24,31,g 


(55638) 


2002 VE 95 


Plutino 


5.143+0.062 


1.790+0.040 


Methanol 


24, g 


(58534) 


Logos 


Cold Classical 


6.759+0.181 


1.653+0.150 




7, 22 


(59358) 


1999 CL 158 


Scattered Disk Object 


6.653+0.090 


1.190+0.072 




8 



References: see AppendixlAl 
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Table A.l cont'd 



Object 

(60454) 2000 CH I05 
(60458) 2000CM 114 
(60558) Echeclus 
(60608) 2000 EE 173 

(60620) 2000 FD 8 

(60621) 2000 FE 8 
(63252) 2001 BL 4 i 
(65489) Ceto 
(66452) 1999 OF 4 
(66652) Borasisi 
(69986) 1998 WW 24 
(69988) 1998 WA 31 
(69990) 1998 WU31 
(73480) 2002 PN 34 
(79360) 1997 CS 29 
(79978) 1999 CCi 58 
(79983) 1999 DF 9 
(80806) 2000 CM 105 
(82075) 2000 YW 134 
(82155) 2001 FZ n3 
(82158) 2001 FP 185 
(83982) Crantor 
(84522) 2002 TC 302 
(84719) 2002 VR I28 
(84922) 2003 VS 2 
(85633) 1998 KR 65 
(86047) 1999 OY 3 
(86177) 1999 RY 215 
(87269) 2000 00 67 
(87555) 2000 QB 243 
(88269) 2001 KF 77 
(90377) Sedna 
(90482) Orcus 
(90568) 2004 GV 9 
(91133) 1998 HK 15 i 
(91205) 1998 US 43 
(91554) 1999 RZ 215 
(95626) 2002 GZ 32 
(118228) 1996 TQ 66 

(118378) 1999 HT„ 

(118379) 1999 HC J2 
(118702) 2000 OM 67 
(119068)2001 KC 77 
(119070) 2001 KP 77 
(119315) 2001 SQ 73 
(119473)2001 U0 18 
(119878)2002 CY 224 
(119951) 2002 KX I4 
(120061)2003 COi 
(120132) 2003FY 128 
(120181)2003 UR 292 
(120216) 2004 EW 95 
(121725) 1999XX I43 
(126619) 2002 CX 154 
(127546) 2002 XU 93 
(129772) 1999 HR„ 
(130391)2000 JG 8i 
(134860) 2000 OJ 67 
(135182) 2001 QT 322 
(136108) Haumea 
(136120) 2003 LG 7 



References: see Appendix|Ap 



Dynamical Class 
Cold Classical 
Scattered Disk Object 
Centaur 

Scattered Disk Object 
Resonant (7:4) 
Resonant (5:2) 
Centaur 

Scattered Disk Object 
Cold Classical 
Cold Classical 
Plutino 

Resonant (5:2) 
Plutino 

Scattered Disk Object 
Cold Classical 
Resonant (12:5) 
Hot Classical 
Cold Classical 
Resonant (8:3) 
Scattered Disk Object 
Scattered Disk Object 
Centaur 

Scattered or Detached KBO 

Plutino 

Plutino 

Cold Classical 
Scattered Disk Object 
Scattered Disk Object 
Scattered Disk Object 
Scattered Disk Object 
Centaur 
Detached KBO 
Scattered Disk Object 
Hot Classical 
Plutino 
Plutino 

Scattered Disk Object 

Centaur 

Plutino 

Resonant (7:4) 

Scattered Disk Object 

Scattered or Detached KBO 

Resonant (5:2) 

Resonant (7:4) 

Centaur 

Plutino 

Resonant (12:5) 
Scattered Disk Object 
Centaur 

Scattered Disk Object 
Scattered Disk Object 
Plutino 
Centaur 

Scattered or Detached KBO 
Scattered Disk Object 
Resonant (7:4) 
Resonant (2:1) 
Cold Classical 
Scattered Disk Object 
Resonant(12:7) 
Resonant (3:1) 



H R (a) 
6.363±0.077 
6.954±0.044 
9.669±0.090 
8.028+0.107 
6.344+0.061 
6.510+0.062 
11.273+0.065 
6.205+0.060 
6.255+0.090 
5.420+0.051 
7.964+0.096 
7.303+0.149 
7.988+0.200 
8.487+0.046 
5.068+0.085 
5.409+0.091 
5.797±0.110 
6.302+0.030 
4.429+0.064 
5.811+0.027 
5.940+0.053 
8.693+0.057 
3.682+0.067 
5.005+0.040 
3.794+0.070 
6.599+0.073 
6.293+0.055 
6.736+0.114 
9.057+0.170 
8.439+0.119 
10.038+0.020 
1.120+0.088 
1.991+0.054 
3.786+0.080 
6.937+0.076 
7.852+0.050 
8.072+0.079 
6.603+0.131 
7.245+0.195 
6.906+0.040 
7.611±0.170 
7.075+0.036 
6.822+0.030 
6.873+0.305 
8.857+0.069 
7.804+0.506 
5.871+0.056 
4.349+0.124 
9.134+0.140 
4.486+0.053 
7.093+0.100 
6.309+0.050 
8.586+0.096 
7.178+0.075 
7.942+0.019 
7.172+0.150 
7.748+0.056 
6.001+0.120 
7.752+0.320 
0.205+0.011 
8.322+0.049 



B-R 
1.699+0.083 
1.240+0.040 
1.376+0.072 
1.164+0.032 
1.806+0.113 
1.230+0.027 
1.199+0.045 
1.420+0.040 
1.830+0.095 
1.610+0.050 
1.235+0.152 
1.412+0.127 
1.225+0.086 
1.280+0.020 
1.746+0.077 
1.566+0.100 
1.630+0.078 
1.980+0.230 
1.417+0.077 
1.418+0.030 
1.402+0.055 
1.864+0.044 
1.741+0.048 
1.540+0.040 
1.520+0.030 
1.727+0.144 
1.055+0.050 
1.151+0.183 
1.702+0.092 
1.088+0.094 
1.810+0.040 
1.874+0.115 
1.042+0.037 
1.470+0.040 
1.240+0.064 
1.185+0.102 
1.346+0.132 
1.199+0.075 
1.881+0.144 
1.830+0.100 
1.384+0.214 
1.290+0.040 
1.470+0.010 
1.720+0.319 
1.130+0.020 
2.079+0.376 
1.680+0.100 
1.660+0.040 
1.240+0.040 
1.650+0.020 
1.690+0.080 
1.080+0.030 
1.734+0.145 
1.470+0.128 
1.200+0.020 
1.450+0.156 
1.417+0.060 
1.720+0.078 
1.240+0.060 
0.973+0.024 
1.271+0.091 



Spectral features 



Water ice 



Water ice 
Featureless 



Methanol 
Water ice 

Water ice 



Methane 
Methane 
Featureless 



Featureless 
Water ice 



Water ice 



References 

28~ 
25 

18, 24 
18, 25 
18, 28 
8, 25 
25,28 
25, g 
28 

16, 35 
8, 28 
28 
28 
25, j 

6, 7, 14, 22, k 
8, 10, 24 
8 

27 

21,25, 28,30,31 

25,28 

25, 30 

25, 28, 33, 34, h 
21,24, 31, g 
24 

24, g 
18, 19 
8, 9, 18 
16, 18 
21,25 
15,28 
25 

21,24, 36,1 
24, 36, m 

24, h 
8, 16 
28 
18 

25, 30, 33 
6,7 

27 
18 
21 
25 

28, 30 
25, 31 
30 
31 

24, 37, h 
25,27 

21,g 

24 

24 

8, 28 

31 

21 

16 

This work 
8 

37 

38, 39, n 
This work 
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Table A.l cont'd 



Object 


Dynamical Class 


H R (a) 


B-R 


Spectral features References 


(136199) Eris 

\ a *j \J a y y j i — > i a lj 


Scattered or Detached KBO 


-1 511+0 033 

A a A A _!_ \J m \y ^y ^y 


1.207+0.088 


IVTpthanp 

IVlVlllilllVw 


24, 36, o 


(136204) 2003 WL 7 


Centaur 


8.670±0.070 


1.230+0.040 




24 


(136472) Makemake 


Hot Classical 


-0.317+0.024 


1.332+0.029 




36,p 


(137294) 1999 RE^is 

\ a *j i x**yf r / i y y y ivj-/^ y 


Cold Classical 


a 091+0 073 

yy - \y y a _!_ yy ■ yj i ^ 


1 700+0 148 

a a / vy vy _!_ vy a a rvy 




18 


(137295) 1999 RB 2 i 6 


Resonant (2: 1) 


7.668+0.096 


1.419+0.142 




18 


(138537) 2000 OK 67 


Cold Classical 


6.093+0.083 


1.540+0.094 




8 


(144897) 2004 UXm 

y a i ruy / / vy vy r vy j. v | yj 


Hot Classical 


4.216+0.087 


1 530+0 020 

i iyy u —i— vy a vy vy 


IVACLlltllUJl 


37 i 


(145480) 2005 TBion 

V -*- * ~y ru vy / vy vy *y x Ay i y{ } 


Detached KBO 


3 949+n 085 


1 540+0 030 

a a *y rvy —i— vy a vy *y vy 




21 


(148209) 2000 CR105 


Detached KBO 


6.191+0.073 


1.273+0.068 




21, 25 


(148780) Altjira 


Hot Classical 


5 885+0 320 

^ ' « yj yj - ' _!_ v./ . ^7 v./ 


1.640+0.170 




30 


(149560) 2003 OZai 

y a r y vy vy j ^ vy vy y v^/_iy | 


Scattered Disk Object 


8 302+0 028 


1 305+0 048 

i vy —i— vy a vy rvy 




Thi*; work" 

x ULLI3 WU1IV 


(168703) 2000 GPi« 


Scattered Disk Object 


5 795+0 061 


1 160+0 057 

a a a vy vy _!_ vy a vy y / 




8 


(181708) 1993 FW 


Hot Classical 


6.572+0.105 


1.625+0.110 




1, 17, 19, 22 


(181855) 1998 WT^ 


Hot Classical 


7.443+0.079 


1.247+0.140 




28, 40 


(181867) 1999 CVns 

y 1 u 1 u vy / / 1 y 7 y v — * y 1 1 


Resonant (7i3)? 


7 067+0 163 

/ ■ \y vy i _!_ yy ■ a yy ^ 


2 130+0 090 

_ a a *y vy _!_ v./ . vy y vy 




27 


(181868) 1999 CGno 

\ 1 u 1 uuu / 1 y y y v — * vy 1 j y 


Scattered Disk Object 


1 004+0 040 

/ . vy v./ r _!_ v / ■ v./ r yj 


1 530+0 080 

iij^v —i— vy a vy vy vy 




27 


(181871) 1999 COist 

\ 1 u 1 u / a / 1 y y y v. 7 l j 


Cold Classical 


g 607+0 030 


1 940+0 090 

a a y rvy —i— vy . vy y vy 




27 


(181874) 1999 HWn 

y a vy a vy / r / 1 y y y x x t t i j 


Scattered or Detached KBO 


g 706+0 062 

vy - / v./ v./ _!_ v./ . v./ yj x— 


1 323+0 043 

a ■ *y — ^ ' _!_ vy a vy r^y 




21, 27 


(182397) 2001 OWw? 

1 A vy ^/ 1 J _V/V/ A V^ T f 2yy / 


Resonant (9:4) 


6 660+0 064 

\y m\y\y\y _!_ V./ . V/ V.' 1 


1 600+0 070 

a a vy vy vy _!_ vy a vy / vy 




21 


(182934) 2002 GJ« 

y 1 u ^ 7 *7 r / ^ vy vy ^ vy y j 


Hot Classical 


5 469+0 187 

~y m r yy y _!_ yy ■ i u / 


1.678+0.261 




30, 31 


1993 RO 


Plutino 


8 492+0 113 


1 385+0 154 

i ■ y u -y _!_ vy a a r 




1, 9 


1994 EV3 


Cold Classical 


7.110+0.072 


1 732+0 167 

a , / ^ ' _ _!_ vy a a vy / 




1 18 27 


1994 TA 


Centaur 


11.421+0.126 


1 930+0 155 




9, 7 


1995 HMs 


Plutino 


7 849+0 109 

/ a \y ry _!_ yy ■ a \y ^/ 


1.010+0.192 




6, 22 


1995 WY 2 


Cold Classical 


6 864+0 110 

yy - yy yy r _!_ yy . a a yy 


1 655+0 278 

i iUyy —i— vy a — / vy 




1 7 


1996 RO™ 


Hot Classical 


ff 903+0 092 

\y • y \J *y _!_ v / . yj y 


1 523+0 156 

a a *y — y _!_ Vi ly u 




7, 10 


1996 RR™ 

A_yyvy xvxv^y 


Plutino 


6.622+0.143 


1 868+0 130 

a a vy vy vy _!_ vy ■ a .y vy 




7, 9, 18 


1996 TK^ 


Cold Classical 


a 190+0 116 

yy m ± y yy _i_ yy . a a yy 


1 666+0 088 

a a vy vy vy _!_ vy a vy vy vy 




7 8 9 

f . O, 7 


1996 TSf,A 


Hot Classical 


5 947+0 130 


1 665+0 157 




6, 7, 12 


1997 CVto 

1 y y > v^ v ^y 


Hot Classical 


7 154+0 030 

/ a a ^y r _!_ vy ■ vy ^y v_/ 


1.860+0.022 




19 


1997 QH 4 


Hot Classical 


a 996+0 136 

vy . y y v.' _i_\y . a ^ vy 


1 731+0 168 

x . / «y x _!_ vy a a vy u 




7 9 10 18 


1997 RT 5 


Hot Classical 


7.117+0.140 


1 549+0 162 

A a *y V _!_ Vy a A V <Cm 




18 


1997 SZm 

I y y 1 wy l—i 1 yj 


Resonant (2: 1) 


8 100+0 104 

\j m a vy vy _!_ vy ■ a vy r 


1 790+0 085 

a a / y vy _!_ vy a vy vy *y 




9 


1998 FSi^/i 


Hot Classical 


6 717+0 105 

\j , I ± I _!_ vy . a vy,y 


1 516+0 057 

i .y i u -!_ vy .yj^j i 




19 22 


1998 HLisi 

Ay y u x aaj 1 j 1 


Hot Classical 


8.120+0.149 


1.190+0.284 




27, 40 


1998 KG« 


Cold Classical 


6.125+0.110 


1 602+0 158 

a a vy vy ^ _!_ vy a a vy 




16 18 


1998 KS« 

17 7 vy Ai.wyfx'^ 


Cold Classical 


7 166+0 040 

/ a a vy vy _!_ vy > vy r vy 


1 730+0 045 

a a / *y vy _!_ vy a vy r*y 




19 


1998 UR/n 

I y y vy vy xViA_^ 


Plutino 


8 083+0 132 

v7 - vy y.> ^ ' _!_ vy ■ a *y — 


1 390+0 113 

i iyy u —i— vy ■ a a 




10 


1998 WS11 

I y y vy y y v. 7 j 


Plutino 


7 952+0 186 

/ ■ y _ _!_ vy ■ i u yy 


1 315+0 075 

a a *y a *y _!_ vy ■ vy / *y 




?8 


1998 WVod 


Cold Classical 


7 126+0 067 

/ a a _ vy _!_ vy ■ vy vy / 


1.270+0.032 




9 


1998 WV 3 | 


Plutino 


7 627+0 069 

/ a w _ / _!_ yy * yy vy v 


1 349+0 096 

a a ^ ' r y _!_ v./ . v./ y vy 




10, 28 


1998 

1 v .7 vy y y z ^-Z^r 


Cold Classical 


6 241+0 099 

vy > •— - r a _!_ v./ « v./ y y 


1 790+0 071 

i ■ / y u —i— vy a vy / a 




9 


1998 WZd, 

1 y y vy tt ^— ' j 1 


Plutino 


8.044+0.102 


1 263+0 089 

a a _ vy -. ' _!_ v./ . vy yj y 




28 


1998 XYos 

1 y y \J 1 v 1 y^ 


Scattered or Detached KBO 


6 438+0 143 

vy > r „ / vy — i— vy ■ a r*y 


1.580+0.212 




14 


1999 CBno 

1 7 7 7 v^jj 1 1 y 


Hot Classical 


6 740+0 050 

vy ■ / rvy _!_ vy ■ vy *y vy 


1 926+0 095 

i . y _ vy _!_ v./ . v / y ^ / 




?8 


1999 CD,« 


Resonant (7:4) 


4.837+0.111 


1 384+0 116 

a a «y vy r_!_ vy a a a vy 




8, 10, 40 


1999 CFhq 


Scattered or Detached KBO 


6 982+0 084 

vy - y vy - _!_ vy . vy vy r 


1.424+0.072 




27, 25 


1999 CJ110 

1 y y y v — ^ y j jy 


Cold Classical 


6 695+0 210 

vy > vy ^ ^y _!_ vy . - a vy 


2.070+0.220 




27 


1999 CMnu 

1 y y y viYi 1 1 y 


Cold Classical 


7 356+0 060 

/ a *y *y vy _!_ vy ■ vy vy vy 


1.780+0.170 




27 


1999 COm 


Hot Classical 


6 682+0 050 

vy a vy vy _ _!_ vy a vy *y vy 


1 350+0 070 

iijjv —i— vy a vy / vy 




27 


1999 CXm 

1 7 7 y v — / 1 j 1 


Resonant (5:3) 


6 914+0 087 

vy a y a r _!_ vy a vy vy / 


1 637+0 118 

i iuJ / —i— vy a a a vy 




?8 


1999 GS,^ 


Hot Classical 


6 230+0 020 

vy a — ^ / vy _!_ vy a vy ^ vy 


1 760+0 070 

a a / vy vy _!_ vy a vy / vy 




27 


1999 HSn 


Cold Classical 


6.344+0.081 


1.845+0.099 




16 19 28 35 


1999 HVn 


Cold Classical 


7 003+0 050 

/ a vy vy *y _!_ vy a vy ^y vy 


1 700+0 063 

Aa / vy vy — i— vy a vy vy y 




19 


1999 JD 132 


Hot Classical 


5.983+0.020 


1.590+0.090 




27 


1999 OE 4 


Cold Classical 


6.887+0.193 


1.832+0.147 




28 


1999 OJ 4 


Cold Classical 


6.899+0.060 


1.675+0.077 




28 


1999 OM 4 


Cold Classical 


7.521±0.100 


1.739+0.170 




18 


1999 RJ 2 i5 


Scattered Disk Object 


7.881+0.103 


1.221±0.175 




18 


1999 RX 214 


Cold Classical 


6.385+0.050 


1.647+0.070 




28 


1999 RY 2I4 


Hot Classical 


7.006+0.040 


1.258+0.085 




28 



References: see Appendix|A*| 
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Table A.l cont'd 



Object 


Dynamical Class 


H R (a) 


B — R Spectral features References 


1999 TR n 


Plutino 


8.063±0.140 


1 770+n 106 


9 


2000 AF 25 , 


Scattered Disk Object 


5.682±0.030 


1 78O+O 060 


27 


2000 CG105 


Hot Classical 


6.469±0.293 


1 170+O 170 

i. . J. / UXU. J. / \J 


27, 40 


2000 CJios 


Hot Classical 


5.687±0.066 


1 760+0 1 06 


31 


2000 CL 104 


Cold Classical 


6.394±0.086 


1 851+0 192 


18 


2000 CLi 05 


Cold Classical 


6.761+0.060 


1 520+0 090 


27 


2000 CN105 


Cold Classical 


5.286±0.160 


1 790+0 1 28 


3 1 

.J 1 


2000 CO 105 


Hot Classical 


5.619+0.124 


1 520+0 180 

_L ■ J-^\J -!—\J .1UU 


27 


2000 CQios 


Scattered Disk Object 


5.996+0.054 


1 1 07+0 043 


25, 28 


2000 FS 53 


Cold Classical 


7.165+0.124 


1 786+0 095 


19 27 

17, Z. / 


2000 FZ 53 


Centaur 


11.103+0.165 


1 170+0 050 


25 


2000 KK4 


Hot Classical 


5.982+0.103 


1 550+O 050 


19 


2000 PE 30 


Scattered Disk Object 


5.867+0.110 


1 1 32+0 084 


15 16 21 


2000 YB 2 


Scattered Disk Object 


6.436+0.084 


1 500+0 1 34 


31 


2001 FM194 


Scattered Disk Object 


7.453+0.159 


1 1 qo+0 040 


25 


2001 HY 65 


Hot Classical 


6.041+0.064 


1 SI 0+0 09? 


3 1 


2001 HZ 58 


Cold Classical 


6.158+0.053 


1 640+0 085 


31 


2001 KA 77 


Hot Classical 


5.050+0.089 


1 812+0 122 

1 . O 1 _ _!_ V / ■ I 


8, 28, 30 


2001 KB77 


Plutino 


7.349+0.078 


1 390+O 1 30 

J. • J* \J -!—\J . M-*J\J 


24 


2001 KD 77 


Plutino 


5.928+0.096 


1 763+0 060 

i ■ / \J*J -!—\J .\J\J\J 


8, 28 


2001 KG 77 


Scattered Disk Object 


8.340+0.120 


1 240+0 070 


25 


2001 KY 76 


Plutino 


6.689+0.380 


1 960+0 291 


30 


2001 QC 298 


Hot Classical 


6.381+0.174 


1 030+0 098 

J. .UJUXU.U7U 


31 


2001 QD 298 


Hot Classical 


6.185±0.170 


1 640+0 1 58 


30 


2001 QF 298 


Plutino 


5.119+0.118 


1 051+0 085 


15, 24, 30 


2001 QR 322 


Neptune Trojan 


7.828+0.010 


1 260+0 036 


41 


2001 QX 322 


Scattered Disk Object 


6.144±0.146 


1 752+O 280 

_L . / JZ.J_U.i-OU 


25, 31 


2001 QY 297 


Cold Classical 


5.151+0.231 


1 561 +0 1 77 


15 30 35 


2001 RZ143 


Cold Classical 


6.241+0.123 


1 <;qo+0 191 

J. .J7uxu. J. y J. 


31 


2001 XZ 255 


Centaur 


10.800+0.080 


1 910+0 070 

1 .71UXU.U / \J 


25 


2002 DH 5 


Centaur 


10.115+0.100 


1 054+0 075 

1 .UJtXU.U / 


28 

Z.O 


2002 GB 32 


Scattered Disk Object 


7.638+0.019 


1 390+O 020 

i. . J 7UXU. V/Z.V/ 


21 


2002 GF 32 


Plutino 


5.973+0.210 


1 7fi5+0 1 34 

1 . / UJXU. 1 Jt 


30 


2002 GH 32 


Hot Classical 


6.098+0.201 


1 509+0 1 60 

1 .JU7XU. 1 \J\J 


30 31 


2002 GP 32 


Resonant (5:2) 


6.580+0.162 


1 386+0 162 

J. . JOUXU. J. \JZ, 


30, 35 


2002 GV 32 


Plutino 


6.886+0.199 


1 860+0 122 

JL . UUUXU. J. Z.Z. 


30 


2002 MS 4 


Resonant (18:11) 


3.333+0.040 


1 070+0 040 


24 


2002 VT 130 


Cold Classical 


5.426+0.092 


2 010+0 233 

Z..U1UXU.Z.JJ 


31 


2002 XV 93 


Plutino 


4.434+0.040 


1 090+0 030 

1 .U7UXU.UJU 


24 


2003 AZ 84 


Plutino 


3.537+0.053 


1 052+0 057 Methanol 


24 31 33 h 


2003 FZ 129 


Scattered or Detached KBO 


6.983+0.038 


1 320+0 040 

J. . JZ/UXU. \J^\J 


21 


2003 HB 57 


Scattered or Detached KBO 


7.389+0.028 


1 310+0 030 

1 .J lUXU.UJU 


21 


2003 QA 92 


Scattered Disk Object 


6.367+0.240 


1 670+0 020 

J. .VJ / \J -1—\J . V/Z.V/ 


37 


2003 QK 91 


Scattered or Detached KBO 


6.966+0.036 


1 370+O 040 

i. .J / UXU.UtU 


21 


2003 QQ 9[ 


Scattered Disk Object 


7.624+0.280 


1 1 80+0 050 


37 

J 1 


2003 QW 90 


Hot Classical 


4.730+0.057 


1 780+0 092 

Jl . / OUXU.U7Z, 


31 


2003 TH 58 


Plutino 


6.940+0.056 


990+0 07 1 

U.77V/XU.U / Jl 


40 


2003 UZ 117 


Hot Classical 


4.920+0.083 


990+O 050 Water ire 

U.7 7UXU.UJV/ VV illL 1 l^C 


24, q 


2003 YL 179 


Cold Classical 


7.482+0.300 


1 260+0 090 


37 


2004 OJ14 


Scattered or Detached KBO 


6.991+0.028 


1 420+0 030 


21 


2004 UP10 


Neptune Trojan 


8.651+0.030 


1 1 60+0 064 


41 


2004 XR I90 


Detached KBO 


3.937+0.036 


1 240+0 040 

1 .Z.~U-LU.U~U 


21 


2005 CB 79 


Hot Classical 


4.375+0.028 


1 090+0 028 Water ice 


40, q 


2005 E0 297 


Resonant (3:1) 


7.221+0.047 


1.320+0.050 


21 


2005 GE 187 


Plutino 


7.192+0.097 


1.740+0.112 


40 


2005 PU 2I 


Scattered Disk Object 


6.091+0.019 


1.790+0.020 


21 


2005 SD 278 


Scattered or Detached KBO 


5.915+0.019 


1.530+0.020 


21 


2005 TN 53 


Neptune Trojan 


9.027+0.040 


1.290+0.106 


41 


2005 TO v4 


Neptune Trojan 


8.426+0.030 


1.340+0.078 


41 


2006 RJ103 


Neptune Trojan 


7.400+0.023 


1.903+0.044 


This work 


2006 SQ 372 


Scattered Disk Object 


7.709+0.049 


1.712+0.093 


21, This work 



References: see , 
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Table A.l cont'd 



Object 


Dynamical Class 


H R (a) 


B — R Spectral features References 


2007 JJ 43 


Hot Classical 


4.044±0.019 


1 610+0 0?0 


21 


2007 JK43 


Scattered Disk Object 


7.028±0.017 


1 400+0 027 


This work 

1 111.) VVWXJV 


2007 NC 7 


Scattered Disk Object 


8.068±0.018 


1.282+0.028 


This work 


2007 RH 283 


Centaur 


8.435+0.039 


1.237+0.069 


This work 


2007 TG422 


Scattered Disk Object 


6.186±0.010 


1.390+0.040 


21 


2007 UM126 


Centaur 


10.161±0.042 


1.080+0.096 Water ice 


This work, i 


2007 VJ305 


Scattered Disk Object 


6.713+0.028 


1.440+0.030 


21 


2008 FC 76 


Centaur 


9.181+0.039 


1.756+0.024 Methanol 


This work, i 


2008 KV 42 


Scattered Disk Object 


8.564+0.056 


1.290+0.060 


21 


2008 OG 19 


Scattered or Detached KBO 


4.612+0.013 


1.470+0.010 


21 


References: see 


Appendix [X| 
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